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ABSTRACT: The present study investigates the creep in a thick-walled composite cylinders made 
up of aluminum/aluminum alloy matrix and reinforced with silicon carbide particles. The distribution 
of SiCp is assumed to be either uniform or decreasing linearly from the inner to the outer radius of 
the cylinder. The creep behavior of the cylinder has been described by threshold stress based creep 
law with a stress exponent of 5. The composite cylinders are subjected to internal pressure which is 
applied gradually and steady state condition of stress is assumed. The creep parameters required to 
be used in creep law, are extracted by conducting regression analysis on the available experimental 
results. The mathematical models have been developed to describe steady state creep in the composite 
cylinder by using von-Mises criterion. Regression analysis is used to obtain the creep parameters 
required in the study. The basic equilibrium equation of the cylinder and other constitutive equations 
have been solved to obtain creep stresses in the cylinder. The effect of varying particle size, particle 
content and temperature on the stresses in the composite cylinder has been analyzed. The study 
revealed that the stress distributions in the cylinder do not vary significantly for various combinations 
of particle size, particle content and operating temperature except for slight variation observed for 
varying particle content. Functionally Graded Materials (FGMs) emerged and led to the development 
of superior heat resistant materials. 

Keywords - FGM(Functionally graded material), Creep law, Von-Mises criterion, Creep Parameters, 



I. Introduction 

There is a development of new class of new material known as composite material, due to 
improvement in performance of structural components in terms of weight reduction and high strength. 
Nowadays composite materials are used very widely because they have desirable properties and are suitable in 
different situations and served their purpose very well. They are made up of two or more constituent which are 
in soluble and mixed at microscopic level. These are heterogeneous materials and generally used to ingredients; 
a binder or matrix and reinforcement (Alman, 2001). The reinforcement is stronger and hard component where 
as matrix act as a binder and helps reinforcement to held in its fixed place. The matrix shields the reinforcement 
which may be breakable like long fibers used with plastics to make fiber glass. 

The composite material has high specific modulus and specific strength than that of monolithic 
materials. Specific modulus is the ratio of Young’s modulus to density where as specific strength is the ratio of 
strength to the density of material for an example, the strength of graphite-epoxy unidirectional composite 
material is same as that of steel but specific strength of this composite material is three time more than steel and 
hence reducing cost of material (Kaw, 1997). 



II. Functionally Graded Materials 

There development of new materials for components exposed to elevated temperature and high 
pressure due to advancement in technology. The conventional metals alone may not survive under high thermo- 
mechanical loading. Thus, a new concept of Functionally Graded Materials (FGM) evolved which led to the 
development of superior heat resistant materials. Such new development materials have capability to resist the 
extreme service conditions. Functionally graded material is a type of multiphase material in which the material 
properties change continuously along the position to meet the design requirements of the components. There 
also achieve the purpose of optimizing the structure (Gupta el al, 2005). Functionally graded materials are 
extensively used in applications such as aerospace, nuclear reactors, internal combustion engines etc. (Noda et 
al, 1998). 

The concept of FGM was first proposed by Japanese researchers in 1984-85 while they working on 
space plane project for advanced material in aerospace application. The body of the space craft will be undergo 
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a very high temperature of about 1700°C with a thermal gradient of 1000°C between inside and outside of the 
space craft. There was no monolithic material that can withstand such extreme temperature conditions. 
Therefore, the researchers proposed a concept to fabricate a material by gradually varying the material 
composition that led to improvement in both thermal resistance and mechanical properties. 



III. Literature Review 

The present study deals with creep behavior in thick-walled functionally graded cylinders exposed to 
operate under high service temperature and pressure. Cylinder has many industrial applications such as pressure 
vessels, emergency breathing cylinders, cylinders for aerospace industries, nuclear reactors and military 
applications, composite rotors (flywheel) etc. In most of these applications such as pressure vessel for storage of 
gases or a media transportation of high -pressurized fluids and composite rotors, the cylinder has to operate 
under severe mechanical and thermal loads that causes creep, and reduced service life. The high creep strength 
along with high specific strength offered by aluminum based metal matrix composites make them a suitable 
choice for these applications. 

3.1 Yield Criteria 

Design of any successful structural components depends upon efficient and safe use of materials. 
Therefore, the ability to predict the yield strength of a material i.e. the onset of plastic yielding, is of great 
practical importance for design purposes. In the field of plasticity, it is important to establish mathematical 
equations for determining the conditions at which yielding starts when a material is under any possible 
combination of loads. Further, in uniaxial loading the plastic flow begins at yield stress as in the case of tensile 
test. The yielding under a combined state of stress can be related to some combination of principal stresses. 

The yield criteria may be expressed in terms of specific quantities, such as the stress state, the strain 
state, strain energy, or others. A yield criterion is usually expressed in mathematical form by means of a yield 

function f (cr , Y) , where CT ; is the state of stress and Y is the yield strength in uniaxial tension (or 



compression). The yield function is defined such that the yield criterion is satisfied when f (CT j] , Y) = 0 . When 

/ (cr , Y) < 0 , the stress state is elastic. The condition /(<T-,F)>0is undefined. To develop a yield 
function, the components of the multi-axial stress are combined into a single quantity known as the effective 
stress ( <J e ). The effective stress is then compared with the yield stress Y, in some appropriate form, to 
determine if yield has occurred. 

For isotropic materials, the yielding would depend on the principal stresses cr l , <J 2 and cr 3 . For 

designing these materials the well-known Von-Mises (1913) yield criterion, also known as distortion-energy 
criterion, is used and is might be expressed as below. 




-cr 2 ) 2 +(a 2 



cr 3 ) 2 +(cr 3 




Where, (j is equivalent stress. 

For ductile materials the shear stress criterion also known as Tresca criterion, may also be used as this 
involved less mathematical complexity. The governing equation of shear stress criterion is given below, 



The graphical comparison of von-Mises criterion and Tresca criterion is given below in Fig. 
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3.2 Variation In Distribution Of Partite Content 

The silicon-carbide particles are assumed to vary linearly from the inner radius a and outer radius b, in 
the thick-walled functionally graded cylinder. As a result, the density of material, conductivity and creep 
parameters will also vary with the radius. The particle content (vol %) of silicon carbide V (r) at any radius r is 
assumed to vary according to the following linear equation, 



V(r) = V -V 1 

v ' max i L max mm J 



b — a 



(3.1) 



Where, V max and V mw are the maximum and minimum particle content respectively at the inner and the 
outer radius. In order to keep average particle content (V avg ) in the cylinder to be constant. One has to select 
either V max or V min . 

Using the law of mixture, the density of FG cylinder is obtained as. 



[100 — V(r f\p m 



100 



-v( r )p d _ ,y> d ~pJ v (n 

/ m 



100 



p(r) = 

(3-2) 

here, p m = 2698.9 km/m 3 and p d = 3210 kg/m 3 are respectively the densities of pure aluminum matrix and SiCp 
(Metals Handbook, 1978; Clyne and Withers, 1993). 

Using the value of V (r) from Eqn. (4.1) into Eqn. (4.2), we obtain 

(p d -pj(r, -y 2 r ) 



P(r) = P m +- 



100 



= X —Y.r 



) 

Where, 



X = p + (jj -pj— and Y = Y ^ Pd Pm) 
f m , d > m 100 100 
The average particle content ( V avg ) in the cylinder can be expressed as, 

b 

J 2nrlV (f) dr 

V = 

avg n (b-a)l 



) 

Where l, denotes the axial length of the composite cylinder. 

Substituting the value of V (r) from Eqn. (4.4) and integrating, we obtain, 

3VL, (l-ff 2 )(l-ff)-V max (l-3ip 2 + 2ip 3 ) 



V. = 

min 



(2-3ip + ip 3 ) 



) 

Where y/ = alb 

3.3 CREEP LAW 



(3.3 



(3.4 



(3.5 



In aluminum based composites, undergoing steady state creep, the effective creep rate ( £ e )is related to 



the effective stress (crj through the well documented threshold stress (a 0 ) based creep law given by (Mishra 
and Pandey, 1990; Park et al, 1990; Mohamed el al, 1992; Pandey and Mishra, 1992; Pandey et al , 1994. 



£., - A 



<?., - CTfi 



exp 



RT 



Where A , n, Q, E, R and T denotes respectively the structure dependent parameter, true stress exponent, true 
activation energy, temperature-dependent Young modules, gas constant and operating temperature. 

The creep law given by Eqn. (3.6) may alternatively be expressed as, 
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s e =\M(r)(a e -a 0 (r)Y 



) 

Where 



M = — 
E 



A exp 



zR 

RT 



and the values of stress exponent n is taken as 5. 



(3.7 



The symbols M(r) and a c (r) are known as creep parameters. 

The creep parameters M (r) and <j 0 (r) given in Eqn. (3.7) are dependent on the type of material and are 
also affected by temperature (T) of application. In a composite, the reinforcement size (P) and the content of 
reinforcement (V) are the primary material variables affecting their creep parameters. In the present study, the 
values of M(r) and a n (r) have been extracted from the uniaxial creep results reported for Al-SiCp by Pandey el 
al (1992). Though, Pandey et al (1992) used a true stress exponent of 8 to describe steady state creep in these 
composites 

In order to extract the values of creep parameters for Al-SiCp, the individual set of creep data reported 
by Pandey et al (1992) have been plotted as £ 1 7 5 versus cron linear scales, Figs. 4.1- 4.3. The values of creep 
parameters M and cr 0 have been obtained from the slope and intercept of these graphs and are reported in Table 

3.1. The £ 1/5 versus a plots corresponding to the observed experimental data points of Al-SiCp composites 
(Pandey et al , 1992) for various combinations of particle size, particle content and temperature, show an 
excellent linearity as evident from Figs. 4. 1-4.3. 

The coefficient of correlation for these plots has been reported as to 0.91 to 0.99 as given in Table 3.1. 



Table 3.1: Creep parameters for Al-SiCp composites (Pandey et al, 1992) 



P 

(um) 


T 

(K) 


V 

(vol%) 


M 

(s 1,s /MPa) 


(MPa) 


Coeff. Of 
correlation 


1.7 


623 


10 


0.00435 


19.83 


0.94 


14.5 


623 


10 


0.00872 


16.50 


0.99 


45.9 


623 


10 


0.00939 


16.29 


0.99 


1.7 


623 


10 


0.00435 


19.83 


0.94 


1.7 


623 


20 


0.00263 


33.02 


0.99 


1.7 


623 


30 


0.00227 


42.56 


0.94 


1.7 


623 


20 


0.00263 


33.02 


0.99 


1.7 


623 


20 


0.00414 


29.79 


0.97 


1.7 


723 


20 


0.00592 


29.18 


0.91 



In order to determine the values of creep parameters M and <x 0 for various combinations of P, V and T, 
reported in Table 3.1, the regression analysis has been performed by using Data-fit software. During the 
regression analysis, P, V and T are taken as independent variables and M and a„ are selected as dependent 
variables. To accomplish the task, the creep parameters given in table have been substituted in the constitutive 
creep model, (4.7) to estimate the strain rates corresponding to the experimental stress values reported by 
Pandey et al (1992) for Al-SiCp composite corresponding to various combinations of material parameters and 
temperature as given in Table 4.1. The developed regression equations are given below. 



Ad (r) = 0.02876 



0.00879 

P 



14.02666 

7\r) 



0.032236 

V(r) 



(3.8 

) 

cr o (r) = -0.084 J P-0.02327Xr) + 1.1853V(r) + 22.2 



(3.9 

) 

Where, P, V(r), T(r), M(r) and a„ respectively the particle size, particle content, temperature, creep 
parameter and threshold stress at any radius (r) of the FG cylinder. 

In a FG cylinder, with particle content and temperature are varying along with radius. Therefore, both 
the creep parameters M(r) and cr 0 (r) will also vary along the radius of the composite cylinder. In the present 
study, the particle size (P) is assumed as 1.7 /-tin while the temperature is assumed to vary along the radius with 
maximum temperature, 450 °C at internal surface and minimum temperature, 350 °C at external surface of FG 
cylinder. Therefore, for a given FG cylinder with known particle and thermal gradients both the creep 
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parameters will be functions of radius. The values of M(r) and a„(r) at any radius (r) could be estimated 
respectively from Eqs. (3.8) and (3.9) by substituting the values of particle size P, particle content V(r) and 
operation temperature, T(r) at the corresponding locations. 



K(r) = 



3.4 Temperature Gradient In The Composite Cylinder 

In order to determine the temperature gradient in the composite cylinder, it is assumed that the internal 
and external surfaces of the cylinder are exposed to a temperature of 450°C and 350 °C respectively. To obtain 
the temperature of any radial distance, T(r), the cylinder is divided into ten elementary cylindrical shells of 
equal thickness, which are, joined together and have varying thermal conductivity. For each cylindrical shell, 
the conductivity is assumed to be constant and determined by law of mixture as given below, 

[100-V(r)] K m +V(r)K d 

(3.10) 

100 

Where K m (= 247 W/mK) is the matrix conductivity (Taya and Arsenault, 1989) and K d (= 100 W/mK) is 
the conductivity of SiCp (Clyne and Withers, 1993). 

Due to the temperature difference between the inner and outer surface of the composite cylinder the heat 
flux per unit length {HA), originated, is given by (Holman, 1992), 

H = 27t(T-T) 

/ In(r 2 / >j)/ K l +In(r 3 / r 2 )/ K 2 + In{r n -r n _ x )l K n 

Where 1) (= 450°C) is the temperature at the inner radius, and T„ (= 350°C) is the temperature at the outer 
radius, and K h K 2 K n are thermal conductivities of different elementary cylindrical shells. 

Knowing 7) and 7'„, the heat flux may be estimated and afterwards the temperature of individual cylindrical 
shells may be estimated from the following equation, 



T =T 

± n + 1 ± n 



H 



2n Kl 



In < r n+JO 



(3.12) 



In the present work four different types of cylinders are assumed, as mentioned below in Table 3.2. The 
conductivity of these cylinders at different radius is estimated from Eqn. (3.10), and is also plotted in Fig. 4.4. 



Table 3.2: Description of cylinders taken in present study 



Cylinder 


Particle Content and other Operating Conditions 


Uniform / Non- FGM(Cl) 


Without particle and thermal gradients [Vmax= Vmin=Vavg = 20 vol%; T = 
Tavg=400 °CJ. Tavg is the average of 77 and To. 


Uniform / Non- FGM (C2) 


With thermal gradient (TG) but without particle gradient (PG) [Vmax= Vmin= 
Vavg = 20 vo/%]. 


FGM (C3) 


With particle gradient (PG) but without thermal gradient (TG) [Vmax= 25 vo/%; 
Vmin= 16 vo/%; Vavg = 20 vo/%; T = Tavg=400 °C], 


FGM (C4) 


With particle and thermal gradients [Vmax= 25 vol%; Vmin= 16 vo/%; Vavg = 
20 vo/%]. 



IV. Mathematical Solution 

Let us consider a thick-walled, circular cylinder made of functionally graded Al-SiCp composite. The 
cylinder is assumed to have inner and outer radii a and b respectively and is subjected to internal pressure p 

only. 

The basic strain rate equations are, 
d/ii 
dr 






(4.1) 



e =tL 

& a — 



(4.2) 
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d/u 

Where Jil = is the radial displacement rate and u is the radial displacement. 

dt 

Where £ r and £ (j are radial and tangential strain rates. 

Eliminating, p from Eqs. (4.13) and (4.14), the basic compatibility equation is obtained as, 

rds a 

' L = £ r - £n 

dr 

The cylinder is subjected to following boundary conditions, 

At r = a, CT r = —p 
At r = b, a =0 



(4.3) 



(4.4) 

(4.5) 



Where a, is the radial stress. 

The equilibrium equation of thick cylinder (Bhatnagar and Arya, 1969) can be written as 
da. 



dr 



(4.6) 



Assuming incompressibility condition (i.e. volume constancy condition) for creep rate, we have 
S r +£g + £ z =0 (4.7) 

The constitutive equations of theory of creep (Bhatnagar and Arya, 1969) are given below, 

£,. 



s a =- 



£_ =■ 



A 

A 

2<7 

A 

2(7, 



[(G + H)a r -Ha d -Ga z ] 


(4.8) 


[( F + H)cr 9 -FcT z -HcT r ] 


(4.9) 


[(F + G)a z -Ga r -Fa 0 \ 


(4.10) 



Where, F, G, H are anisotropic constants of the material and £ e is the effective creep rate invariant. 
The von Mises (Dieter, 1998) effective stress invariant <j e is given as below. 



I [ H A - ' ? + G A -<t z ) 2 +F - a of 



1/2 



(4.11) 



The plane strain condition is assumed i.e. strain rate along axial direction is zero in the present 
investigation, 

£ z =0 (4.12) 

From Eqs. (4.13), (4.14), (4.19), displacement rate is, 

jl = - (4.13) 



Where ‘c’ is the constant of integration. 

Substituting Eqn. (4.25) in Eqn. (4.13), we get, 
c 

By substituting Eqn. (4.25) in Eqn. (4.14), we obtain, 
c 

£ e= 2 

r 

Also from Eqs. (4.22) and (4.24), we have 
Ga. + Fa n 



F + G 

By using Eqn. (4.28) in Eqn. (4.23), the effective stress invariant (cr e ) is given by. 



(4.14) 



(4.15) 



(4.16) 
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1 



FG + GH + HF 



1/2 



cr =—^=(cr -cr g ) - — ' (4.17) 

e y/2 l ( F + G ) 

Substituting the values of ( £ r ) and (<j-) from Eqs. (4.26) and (4.28) in Eqn. (4.20), we obtain. 



(a-aff = 



F + G 



FG + HF + GH 



2<X c 



£„ r 



(4.18) 



Substituting the values of (cr e ) and ( £ e ) from Eqs. (4.29) and (4.7) in Eqn. (4.30), we obtain. 



(cj r -(J e )= — 



(4.19) 



Where, 
T(r) = 



2 (F + G) 


n+ 1 
2 n 


~-C Un ~ 


(FG + GH + HF 




M(r) 



P(r) = 



2 (F + G) 






(FG + GH + HF 
Solving Eqs. (4.18) and (4.31), we have. 



cla r _ T(r ) 
dr r 2,n 



P(r) 



da, T(r) P(r) 



dr 



n+2 



By integrating above expression from limit a to r, we get 
T(r) , r P(r) 



C' i \ r) , , f' t'K.r) , , ^ 

= 1 ~Frir- dr+ I -dr + C, 

•> Cl J a y 



(4.20) 



Where C/ is another constant of integration that can be evaluated by applying boundary conditions given by 
Eqn (4.4) in Eqn (4.20) 

Therefore, Eqn (4.20) becomes, 

— • “r-P (4 - 2l » 

Ja Ja 

r n 

By applying another boundary condition given by Eqn. (4.5) in Eqn. (4.20), the constant ‘C’ appearing in 
expression of T(r) can also be evaluated as below, 

— tl/2 



P~ 



c = 



2 (F+G) 



(FG + GH+HF) 



rb CT. 



O) 



.dr 



F + G 



(FG + GH + HF) 



n + 1 

2 n r b 



■j: 



«+2 

r^E _m 



.dr 



0-1 



Using, Eqs. (4.21) into Eqs. (4.19), tangential stress can be evaluated as below, 
TO) 






PP. dr - ™ 



P -P 
Gr) r 



n+2 • J„ • 2 

Ja y 

y. n y.n 

Putting, Eqs. (4.23) and (4.21) in Eqs. (4.16), axial stress can be evaluated 



['T(n , . r P(r) , T(r) P(r ) 



-P 



2 r n 



(4.22) 



(4.23) 



(4.24) 
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4.1 Calculation Procedure 

Based on analysis carried out in previous section, a mathematical model has been developed to 
calculate the steady state creep behavior of the rotating FGM cylinder consists of linearly decreasing SiCp 
particles and subjected to positive thermal gradient. For the purpose of numerical computation, the inner and 
outer radii of the cylinder are taken 25.4 mm and 50.8 mm respectively. The composite cylinder is subjected to 
internal pressure of 85.25 MPa and zero external pressures. The particle size of SiC is taken as U/jm. The 
dimensions of cylinder selected in this study are similar to those used in earlier experimental work (Johnson et 
al, 1961) on thick- walled cylinder made of aluminum alloy. The radial, tangential and axial stresses at different 
radial locations of the cylinder are calculated respectively from Eqs. (4.21), (4.23) and (4.24). The creep 

parameters M(r) and <7 0 (r) required during the computation process are estimated from Eqs. (3.8) and (3.9), for 

the desired combination of particle size, particle content and operating temperature. Further for isotropic 
cylinders, the anisotropic constants becomes equal and their value is taken as unity i.e. F=G=H= 1. In order to 
determine the thermal gradient the temperature at inner surface is taken as 450°C and at outer surface 350°C 
respectively. The thermal gradient required in Eqs. (3.8) and (3.9) is estimated from Eqn. (3.12). 




particle size of SiC Volume fraction of SiC 




operating temperature. 
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V. Results and Discussions 

Based on the analysis, a mathematical model has been developed to obtain the stress distribution in 
functionally graded cylinder which has been obtained for various combinations of particle size, P, particle 
content, V, and temperature, T. The values of creep parameter have been taken from Table 3.1. The results have 
been obtained for uniform composite cylinder and functionally graded cylinders with and without thermal 
gradient, as mentioned in Table 3.2. 

To accomplish this task, the creep rates are estimated from the present analysis for an isotropic 
composite cylinder made of copper for which the results are reported in the literature (Johnson el al, 1961). The 
dimensions of the cylinder, operating pressure and temperature, and the values of creep parameters used for the 
purpose of validation are summarized in Table 5.1. 

Table 5.1: Data used for validation (Johnson et al, 1961) 

Cylinder Material : Copper 

Cylinder dimensions: a = 25.4 mm, b =50.8 mm. 

Internal Pressure = 23.25 MPa, External Pressure = 0 
Operating Temperature = 250 °C 

Creep parameters estimated: M = 3.271xl0' 4 s' I/s /MPa, a 0 = 1 1.32 MPa 



To estimate the values of parameters M and <7 g for copper cylinder, firstly the values of effective stress, 
<T e , have been calculated at the inner and the outer radii of cylinder by substituting the values of (J r ,(J g and 
<7. into von-Mises yield criterion given by Eqn. 4.11, as reported by Johnson et al (1961), at these locations. 
The effective stresse, thus estimated at the inner radius (CT, = 1 89.83 A// 3 !/ and £ e = 2.168x10 8 S ') and 
the outer ( (J e = 1 16 MPd and £ e = 1.128x10 9 S ' ) of the copper cylinder are substituted in creep law, 

Eqn. (3.7) to obtain the creep parameters M and <J 0 for copper cylinder as given in Table 5.1. The material is 
isotropic, so the values of anisotropic constants, F = G = H = 1. 

5.1 Creep Stresses and Creep Rates 

To study the effect of thermal gradient and particle gradient on creep behavior of composite cylinder, 
the steady state stresses, have been obtained in four different composite cylinders as mentioned in Table 4.2, 
and are also plotted in Figs. (5. 1 )-(5.4). The cylinder taken for investigation in the present work (Table 5.1) is 
having same dimensions that of cylinder used by Johnson et al ( 1961 ) in their work. 

The radial stress, Fig. 5.1 remains compressive (negative) throughout the cylinder, with maximum value 
at the inner radius (82.25 MPa) and zero at the outer radius, due to imposed boundary conditions given in Eqs. 
(4.4) and (4.5). The magnitude of radial stress increases throughout the cylinder when thermal gradient is 
imposed and due to this the variation of radial stress becomes almost linear. However, the presence of particle 
gradient leads to reduction in the magnitude of radial stress over the entire radial distance. The simultaneous 
presence of particle and thermal gradient has averaging effect on the radial stress. It can be seen from Fig. 5. 1 
that the maximum variation observed in the magnitude of radial stress is 5.8 MPa at a radius of 35.6 mm of 
composite cylinder. The tangential stress shown in Fig. 5.2 remains tensile throughout and is observed to 
increase linearly with increasing radial distance for uniform cylinder Cl. The imposition of thermal gradient, 
with temperature at inner radius being higher than at the outer radius, in uniform cylinder (C2) leads to decrease 
in tangential stress by 29.2 MPa, near the inner radius but increase in tangential stress by 45 MPa near the outer 
radius as compared to those observed in cylinder Cl. Further, similar to radial stress, the variation of tangential 
stress also becomes linear and steeper with the imposition of thermal gradient alone. By incorporating more 
amounts of silicon carbide particles near the inner radius, as in FG cylinder C3, the tangential stress 
significantly increases by 30.5 MPa near the inner radius but decreases by 47. 1 MPa at the outer radius in 
comparison to uniform cylinder Cl. Further, the tangential stress distribution becomes more uniform thereby 
reducing the in -homogeneity. The region having relatively more particle offers higher stress compared to the 
region containing lesser amount of particle. The similar observation had also been reported by Pandey et al, 
(1992) in their experimental study of Al-SiC composites. When the cylinder (C4) is subjected to particle and 
thermal gradients together, the tangential stress is having averaging effect and its value is in-between the values 
obtained for cylinder C2 and C3. The axial stress in cylinder Cl changes its nature from compressive (negative) 
at the inner radius and it become tensile (positive) at the outer radius as shown in Fig. 5.3. By imposing 
temperature gradient in uniform cylinder (i.e. cylinder C2), the axial stress near the inner radius becomes more 
compressive. Further, the tensile value of axial stress increases near the outer radius compared to that observed 
in uniform composite cylinder C 1 at constant temperature. The imposition of particle gradient alone in uniform 
cylinder (i.e. cylinder C3) decreases the magnitude of compressive value of axial stress near the inner radius 
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when compared with uniform cylinder C 1 . The simultaneous presence of both the gradients superimposes each 
other to result effect similar to that observed for tangential stress in Fig. 5.2. The distribution of von-Mises 
effective stress is plotted in Fig. 5.4. The effect of thermal and particle gradients alone or together on variation 
of effective stress. Fig. 5.4, is almost opposite in nature to those observed for tangential stress in Fig. 5.2. 
Further, the magnitude of effective stress is also increases for the entire radial distance of the composite 
cylinder as compared to tangential stress values. 



VI. Conclusions 

1 ) The radial, tangential and axial stresses increases as we move from inner towards outer radii of cylinder, 
reaches maximum before decreasing near the outer radius. 

2) The stress distribution in the cylinder does not vary significantly for various combinations of particle size 
and operating temperature. The increase in particle content shows sizable variation on stress distribution 
cylinder 

3) The effect of imposing thermal and particle gradients, either alone or simultaneously, on the axial stress is 
almost similar to that observed for tangential stress in the cylinder. 

4) The effective stress is significantly affected by incorporating particle as well as temperature gradient alone 
or simultaneously. The trend of variation in effective stress is opposite to that of tangential stress. 

Figures: - 





Fig 5.1 : Variation of radial stress in composite cylinder. Fig 5.2: Variation of tangential stress in 

composite cylinder. 





Fig 5.3 : Distribution of axial stress in composite cylinder. Fig 5.4: Distribution of von-Mises effective stress. 
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